US008675922B1 


d2) United States Patent 

Watson et al. 


(io) Patent No.: US 8,675,922 B1 

(45) Date of Patent: Mar. 18, 2014 


(54) VISIBLE MOTION BLUR 

(75) Inventors: Andrew B. Watson, Los Gatos, CA 
(US); Albert J. Ahumada, Mountain 
View, CA (US) 

(73) Assignee: The United States of America as 

Represented by the Administrator of 
the National Aeronautics & Space 
Administration (NASA), Washington, 
DC (US) 

( * ) Notice: Subject to any disclaimer, the term of this 

patent is extended or adjusted under 35 
U.S.C. 154(b) by 142 days. 

(21) Appl.No.: 13/444,777 

(22) Filed: Apr. 11, 2012 


Related U.S. Application Data 


(60) 

Provisional application No. 61/520,357, filed on May 


24, 2011. 


(51) 

Int. Cl. 



G06K 9/00 

(2006.01) 


H04N 5/228 

(2006.01) 

(52) 

U.S. Cl. 



USPC 

382/107; 382/255; 348/208.4 

(58) 

Field of Classification Search 


USPC 

382/100, 103, 106-107, 162, 167, 168, 


382/173, 181, 191, 199, 232, 254-255, 
382/274-275, 276, 305, 312; 348/222.1, 
348/208.4; 250/221; 345/102, 690 
See application file for complete search history. 

(56) References Cited 

U.S. PATENT DOCUMENTS 

7,317,445 B2 * 1/2008 Hekstraetal 345/102 

7,783,130 B2 8/2010 Watson 


2005/0168492 Al * 8/2005 Hekstraetal 345/690 

2008/0061220 Al * 3/2008 Machida 250/221 

2008/0170124 Al * 7/2008 Hatanakaetal 348/208.4 

2010/0066850 Al * 3/2010 Wilson etal 348/222.1 


OTHER PUBLICATIONS 

Video Electronics Standards Association (VESA), "Flat Panel Dis- 
play Measurements Standard (FPDM), Version 2.0,” Video Electron- 
ics Standards Association, Milpitas, CA, 2001. 

G. E. Legge and J. M. Foley, “Contrast masking in human vision,” J. 
Opt. Soc. Am., 70, 1458-71, 1980. 

J. R. Hamerly and C. A. Dvorak, “Detection and discrimination of 
blur in edges and lines,” J. Opt. Soc. Am., 71, 448, 1981. 

E. Peli, “Contrast in complex images,” J. Opt. Soc. Am. A, 7, 2032- 
40, 1990. 

A. B. Watson and J. A. Solomon, “Model of visual contrast gain 
control and pattern masking,” J. Opt. Soc. Am. A, 14,2379-91, 1997. 
G. Westheimer, S. Brincat, and C. Wehrhahn, “Contrast dependency 
of foveal spatial functions: orientation, vernier, separation, blur and 
displacement discrimination and the tilt and Poggendorff illusions,” 
Vision Research, 39, 1631-39, 1999. 

(Continued) 

Primary Examiner — Seyed Azarian 

(74) Attorney, Agent, or Firm — Christopher J. Menke; 

Robert M. Padilla 

(57) ABSTRACT 

A method of measuring motion blur is disclosed comprising 
obtaining a moving edge temporal profile r x (k) of an image of 
a high-contrast moving edge, calculating the masked local 
contrast m^k) for r x (k) and the masked local contrast m 2 (k) 
for an ideal step edge waveform r 2 (k) with the same amplitude 
as r ; (k), and calculating the measure or motion blur V as a 
difference function, 

The masked local contrasts are calculated using a set of con- 
volution kernels scaled to simulate the performance of the 
human visual system, and T 1 is measured in units of just- 
noticeable differences. 
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VISIBLE MOTION BLUR 

RELATED APPLICATION 

This application claims priority to U.S. Provisional Patent 5 
Application No. 61/520,357 filed May 24, 2011, which is 
hereby incorporated by reference. 

ORIGIN OF THE INVENTION 

The invention described herein was made by an employee 
of the United States Government and may be manufactured 
and used by or for the Government of the United States of 
America for governmental purposes without the payment of 
any royalties thereon or therefor. 

FIELD OF THE INVENTION 

One or more embodiments of the present invention relate to 
methods for measuring motion blur in imaging systems. 

BACKGROUND 

Motion blur is a significant defect of most current display 
technologies. Motion blur arises when the display presents 
individual frames that persist for significant fractions of a 
frame duration. When the eye smoothly tracks a moving 
image, the image is smeared across the retina during the 
frame duration. Although motion blur may be manifest in any 
moving image, one widely used test pattern is a moving edge. 
This pattern gives rise to measurements of what is called 
moving-edge blur. 

A number of methods have been developed to measure 
moving edge blur, among them pursuit cameras, so-called 
digital pursuit cameras, and calculations starting from the 
step response of the display. These methods generally yield a 
waveform — the moving edge temporal profile (METP) — that 
describes the cross-sectional profile of the blur [1], 

Several methods have also been developed to convert this 
waveform to a single-number metric of motion blur. 
Examples are the Blur Edge Time (BET), Gaussian Edge 
Time (GET), and Perceptual Blur Edge Time (PBET) [1], 
However, none of these metrics attempts to provide a percep- 
tual measure of the amount of motion blur. 

First, none of these metrics takes into account the contrast 
of the edge, and its effect upon perceived blur. In general, blur 
becomes less visible when contrast decreases [2, 3], and the 
apparent width of motion blur declines with reduced contrast 
[4] . Second, contrast of the edge will mask the visibility of the 
blur [5, 6]. Thus a model of blur visibility must take into 
account this masking effect. 

The need to incorporate contrast is especially pressing 
because measurements of motion blur are often made at sev- 
eral contrasts (gray-to-gray transitions) [7, 8]. Those separate 
measurements must then be combined in some perceptually 
relevant way. 

Finally, none of the existing metrics take into account the 
visual resolution of the display (pixel s/degree of visual 
angle). For a given speed in pixels/frame, a higher visual 
resolution will yield a less visible artifact. 

SUMMARY OF THE INVENTION 

A method of measuring motion blur is disclosed compris- 
ing obtaining a moving edge temporal profile r^k) of an 
image of a high-contrast moving edge, calculating the 
masked local contrast in , (k ) for r^k) and the masked local 


2 

contrast m 2 (k) for an ideal step edge waveform r 2 (k) with the 
same amplitude as r^k), and calculating the measure or 
motion blur T* as a difference function, 

'P=5(Ai2j.lm , (k)-m 2 (k) ! 1 U 1,|i . 

The masked local contrasts are calculated using a set of con- 
volution kernels scaled to simulate the performance of the 
human visual system, and 'P is measured in units of just- 
noticeable differences. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 shows an example of a moving edge temporal pro- 
file (METP) for a blurred edge. 

FIG. 2 shows a fit of a cumulative Gaussian curve to the 
waveform of FIG. 1 . 

FIGS. 3 A-C show examples of the center, surround, and 
masking kernels. 

FIG. 4 shows the results of the convolutions of the center 
and surround kernels of FIG. 3 with the waveform of FIG. 2. 

FIG. 5 shows the contrast waveform, local contrast energy, 
and masked local contrast for the waveform of FIG. 2. 

FIG. 6 shows an ideal step edge overlaid on the METP 
waveform. 

FIG. 7 shows the masked local contrast waveforms for the 
two waveforms of FIG. 6. 

FIG. 8 shows the difference between the two masked local 
contrast waveforms of FIG. 7. 

FIG. 9 shows the value of visual motion blur as a function 
of the offset of the ideal step edge waveform from the METP 
waveform. 

DETAILED DESCRIPTION 

Before the present invention is described in detail, it is to be 
understood that unless otherwise indicated this invention is 
not limited to specific imaging systems. 

It must be noted that as used herein and in the claims, the 
singular forms “a,” “and” and “the” include plural referents 
unless the context clearly dictates otherwise. Thus, for 
example, reference to “an imaging system” includes two or 
more imaging systems, and so forth. 

Where a range of values is provided, it is understood that 
each intervening value, to the tenth of the unit of the lower 
limit unless the context clearly dictates otherwise, between 
the upper and lower limit of that range, and any other stated or 
intervening value in that stated range, is encompassed within 
the invention. The upper and lower limits of these smaller 
ranges may independently be included in the smaller ranges, 
and are also encompassed within the invention, subject to any 
specifically excluded limit in the stated range. Where the 
stated range includes one or both of the limits, ranges exclud- 
ing either or both of those included limits are also included in 
the invention. Where the modifier “about” is used, variations 
of ±1 0% are considered to be within the scope of the disclosed 
limit. 
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As used herein, the following symbol definitions apply: 


TABLE 1 


Symbol 

Definition 

Example 

Unit 

k 

integer sample index 

N, N, 

2 2 

dimensionless 

r(k) 

an arbitrary luminance waveform 


relative luminance 

ri(k) 

r(k) for a moving edge (METP) 


relative luminance 

r 2 (k) 

r(k) for an ideal step edge 


relative luminance 

R o 

lower luminance level for step 

50 

relative luminance 

Ri 

upper luminance level for step 

330 

relative luminance 

At 

time between samples 

0.02867 

frames 

Ax 

distance between samples 

0.007167 

degrees of visual angle 

P 

speed of moving edge 

16 

pixels/frame 

V 

visual resolution 

64 

pixels/degree 

P 

center of cumulative Gaussian 


degrees of visual angle 

a 

standard deviation of cumulative Gaussian 

0.0468 

degrees of visual angle 

g 

cumulative Gaussian 


relative luminance 

N 

number of standard deviations for trim 

32 

dimensionless 

N , 

number of samples in waveform after trim 


dimensionless 

h,(k) 

center kernel 


dimensionless 

s c 

scale of center kernel 

2.77 

degrees of visual angle 

h s (k) 

surround kernel 


dimensionless 

S s 

scale of surround kernel 

21.6 

degrees of visual angle 

hjk) 

masking kernel 


dimensionless 

s w 

scale of masking kernel 

10 

degrees of visual angle 

K 

adaptation weight 

0.772 

dimensionless 

R 

mean relative luminance, 
typically (Rq + R x )/2 

190 

relative luminance 

T 

masking threshold 

0.3 

contrast 

S 

sensitivity 

217.6 

dimensionless 

P 

pooling exponent 

2 

dimensionless 


The present invention discloses a perceptual metric for 
motion blur called Visible Motion Blur (VMB). VMB incor- 
porates the three effects described in the Background section: 
contrast, masking, and visual resolution. It is based on the 
Spatial Standard Observer [9]. VMB converts a moving edge 
temporal profile (METP) waveform into a measure of motion 
blur quantified in units of just noticeable differences (JNDs). 
JND is a standard perceptual measure in which one JND is the 
least quantity that can be seen with specified reliability. 

The starting point for the VMB metric is the METP, a 
discrete sequence of relative luminances, which we write here 
as r^k), where k represents an integer sample index, and the 
time between samples is At in units of frames. This waveform 
is a standard physical measurement of motion blur and can be 
acquired in several ways [ 1 ] . These generally involve estimat- 
ing the width of an edge subjected to motion blur. 

The edge can be captured in any of three ways. The first 
method employs a pursuit camera that tracks a vertical edge 
(between two gray levels) as it moves horizontally across the 
screen. The camera is simulating the eye as it pursues the 
moving edge. The result, after averaging over time, is a pic- 
ture of the blurred edge. After averaging over the vertical 
dimension (orthogonal to the motion), a one-dimensional 
waveform representing the cross-section of the blurred edge 
can be obtained. It describes relative luminance (a linear 
function of luminance) as a function of horizontal position in 
pixels. When recorded at several speeds of edge motion, the 
waveforms are usually found to correspond when the hori- 
zontal scale is divided by the speed. Therefore, it is conven- 
tional to rescale the horizontal axis of the profile (pixels) by 
dividing by the speed (pixels/frame) to obtain a wavefonn 
that is a function of time (frames) — the METP. It is also 
conventional to characterize the width of the METP in terms 
of the time interval between 1 0% and 90% points of the curve. 
This quantity is called the blur edge time (BET). 


The second method employs a stationary high-speed cam- 
era. With a sufficiently high frame rate, it is possible to cap- 
ture a sequence of frames, that, with appropriate shifting and 
adding, can also yield a record of the METP. The high-speed 
camera avoids the mechanical challenges of the pursuit cam- 
era. This second method can be called “digital pursuit.” 

The third method employs a fixed non-imaging detector 
such as a photodiode, which measures the luminance over 
time as the display is switched from one gray level to another. 
This temporal step response is then convolved with a pulse of 
duration equal to the hold time (for an LCD, typically one 
frame), to obtain another version of the METP. This third 
method can be called the “temporal step” method. The tem- 
poral step method relies on an assumption that all pixels are 
independent. It has been demonstrated to be accurate in many 
cases, but may fail when motion-dependent processing is 
present. 

An example of an METP is shown in FIG. 1. In this 
example At=0.02867 (i.e., lAs frame). At can be selected so 
that there are at least 1 0 samples across the step in luminance 
so that the blur is well resolved. The data from FIG. 1 will be 
used throughout the exemplary embodiment of the invention 
below. Note that FIG. 1 has a non-zero black-level. This is 
typical of situations where the METP is recorded in a dark 
environment, but the visibility of motion blur is to be esti- 
mated for a lit environment. A suitable “veiling luminance” 
can be added to the METP to accommodate this background 
level. 

.An exemplary embodiment of the calculation of the VMB 
metric is as follows. The first step is to determine the distance 
between samples Ax in units of degree of visual angle. This is 
given by 
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pht 
Ax= — 

V 


(1) 


where p is the assumed speed of edge motion in pixels/frame 
and v is the visual resolution of the display in pixels/degree. 
For example, if p=16 pixels/frame and v=64 pixels/degree, 
then Ax=0.007167 degrees. 

The waveform r(k) consists of a transition between two 
relative luminance levels R 0 and (FIG. 1). (R 0 is non-zero 
in the example of FIG. 1 to show inclusion of veiling lumi- 
nance.) It is useful (although not necessary) to trim the 
sequence r : (k) to the neighborhood of the transition to reduce 
subsequent computations. The waveform can be fitted to a 
cumulative Gaussian 


g(k',p, cr, R 0 , R l ) = R 0 + 



( 2 ) 


20 


where p is the center of the Gaussian, and ct is the width 
(standard deviation). The wavefonn can then be trimmed to 
values of kthat are within Nstandard deviations of the mean, 
that is, a portion of r t (k) is selected for which 25 

\kAx-p\sNo. (3) 

FIG. 2 shows the fit to the example METP wavefonn from 
FIG. 1. The waveform has also been trimmed using N=32, 
and the horizontal coordinate has been converted to degrees 30 
relative to p (i.e., the horizontal axis is kAx rather than k). It 
can also be convenient to make the length of the sequence an 
even number by deleting the last point if necessary. The 
length of the trimmed sequence is N. The waveform shown in 
FIG. 2, plotted against a spatial distance (degrees of visual 35 
angle) instead of time (frames) is sometimes referred to as a 
“Moving Edge Spatial Waveform” (MESP). When this termi- 
nology is used, an METP is converted to an MESP using the 
relationship between Ax and At given in equation 1 . 

The next step is to create three convolution kernels, h c (k), 40 
h,(k), and h„,(k). These are discrete sequences obtained by 
evaluating kernel functions at a discrete set of points with 
k-values matching those of the trimmed METP wavefonn: 

45 

N, N, (4) 

£ = -— I ... -1-1. 

2 2 

Exemplary embodiments of the convolution kernels are 

, 50 

given by 


h c (k) = — sech7r , 

S c s c 



(5) 

55 

( 6 ) 

(7) 

60 


These are called the “center” kernel, the “surround” kernel, 
and the “masking” kernel respectively. These kernels have 
“scales” (i.e., widths in the k-direction) of s c , s s , and s„„ 
respectively, measured in degrees of visual angle. Each kernel 65 
is normalized to have an integral of 1 . The first two can be 
thought of as simulating the processing of the luminance 


6 

waveform by retinal ganglion cells with antagonistic center 
and surround components. Values of about 2.77 min and 21.6 
min (i.e., 2.77/60 and 21.6/60 degrees of visual angle) are 
found to approximate human visual sensitivity. The center 
component incorporates the blur due to the visual optics, and 
possibly further early neural pooling, while the surround 
computes an average of the local luminance, and uses it to 
convert luminance to local contrast [10]. With the range 
defined by Equation 4, the kernels are each of the same length 
as the trimmed sequence. To reduce computation, they can 
alternatively be made of shorter and different lengths, each 
approximately four times its respective scale. Examples of 
the three kernels are shown in FIGS. 3 A-C using a horizontal 
axis scale corresponding to the center one third of that of FIG. 
2 . 

One of ordinary skill will recognize that the functional 
form of all of the kernels can vary provided that they generally 
have the indicated scales and are suitably normalized to have 
an integral of 1 . The surround and masking kernel examples 
use Gaussian waveforms, while the center kernel example 
uses a hyperbolic secant. These produce similar peaked wave- 
forms with differing tail shapes: a Gaussian tail decreases as 
exp(-[kAx] 2 ), while the hyperbolic secant tail decreases as 
exp (-k0 Ax). Other similar peaked waveforms can also be 
used with similar results. A Cauchy or Lorentz waveform has 
tails which decrease as (kA0 x)~ 2 . Similar functional forms 
can be readily devised which decrease as any even power of k 
0 Ax. The special case of “zero” power is also possible using 
a rectangular waveform with a width equal to one over the 
height. The example waveforms given in equations 5-7 are 
generally found to provide good correlation with the charac- 
teristics of the human visual system. 

The trimmed waveform is convolved with the center and 
surround kernels h c and h 5 to yield h c .*r 1 and h^*^, where * is 
the convolution operator. For example. 


kJkFrJk^MKtt-i) Ax. (8) 

In principle, the sum is over all i-values from -oo to oo; in 
practice, it is sufficient to sum over i -values where h c (i) differs 
significantly from zero. FIG. 4 shows the results of these two 
convolutions for the waveform of FIG. 2 and the kernels of 
FIG. 3. The convolution with the center kernel is shown as a 
solid line, and the convolution with the surround kernel is 
shown as a dashed line. 

Next, the local contrast waveform c(k) is computed. c(k) is 
defined by 


c{k) = - 


h c *r\ 
K(h s *ri) + (1 - k)R 


(9) 


where k is an “adaptation weight” parameter and R is the 
mean relative luminance, typically computed as the average 
of the maximum and minimum relative luminances R 0 and 
Rj, as estimated from the fit of the cumulative Gaussian of 
Equation 2. The effective local contrast energy e(k) is com- 
puted using the masking kernel h rn and a masking threshold 
parameter T: 


{ c(k) i 2 

e(k) = A m (A)tj^ — J . 


(10) 
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TABLE 2 


7 

The masked local contrast m(k) is computed as 


m(k) = 


c(k) 

V 1 + e(k) 


(ID 

5 


FIG. 5 shows Cj(k) (solid line), e^k) (short-dashed line), and 
nqfk) (long-dashed line) for the waveform of FIG. 2 and the 
kernels of FIG. 3. This model of masking is similar to that 10 
developed by Ahumada in a study of symbol discrimination 
[11, 12, 13]. In these calculations, the local contrast energy 
e(k) is a measure of the visually effective pattern ensemble in 
the neighborhood of a point k, and it determines the amount of 
masking of nearby contrast patterns. Patterns are less visible 15 
when they are superimposed on other patterns. 

To compute the visibility of the motion blur, we compare 
the masked contrast waveform mj for the test edge as com- 
puted above to the masked contrast waveform m 2 for an ideal 
edge of the same starting and ending luminance. The visible -, 0 
motion blur (VMB) can now be calculated in units of just 
noticeable differences (JNDs) by 

V=S(bxZ. k \m l (k)-m 2 (k)\V) l! V, (12) 

where S and (3 are parameters (a “sensitivity” and a “pooling 
exponent”). The location of the ideal edge is adjusted to find 25 
the minimum value of V [12, 13]. The value of *P will still 
depend on the alignment of the blurred and ideal edges. The 
effective visible difference corresponds to the minimum of *P. 
This can be determined by computing V for various shifts of 
the ideal edge, as described below. 30 

In greater detail, an ideal step edge waveform (defined for 
the k-range given by Equation 4) is given by 

r 2 (*)=2? 0 +(i? 1 -J? 0 )step(*-S), (13) 

where step is the unit step function, and 6 is between 1 and N,. 35 
8 can, but need not, be an integer. The contrast waveforms 
c 2 (k), e 2 (k), and m 2 (k) are computed as above substituting 
r 2 (k) for r^k), and then JND is computed using Equation 12. 
This process is repeated for each possible value of 8 and the 
smallest value of iris selected as the final value of VMB. 40 

FIG. 6 shows the two input waveforms: the example wave- 
form rj(k) (dotted line) and the ideal edge r,(k) (solid line), 
and FIG. 7 shows the masked local contrast waveforms m t (k) 
(dotted line) and m 2 (k) (solid line). FIG. 8 shows the differ- 
ence between the two masked local contrasts. FIG. 9 shows i]i 45 
as a function of the shift 8. The minimum is 7.5 JNDs. In this 
example, the motion blur is calculated to be clearly visible, 
because the VMB is substantially larger than 1 JND. 

There are several adjustable parameters in the calculation 
described above. These parameters can be chosen to mimic 50 
the response of the human visual system to the motion blur of 
a particular test image. The parameters are summarized in 
Table 2. Values that have been tuned to provide a good match 
to human perception are provided as “example values.” How- 
ever, the parameters can vary over the “usable range” to 55 
provide alternate embodiments of the present invention. The 
calculation of VMB incorporates several important features 
of human contrast detection: light adaptation (in the conver- 
sion to contrast), a contrast sensitivity function (via convolu- 
tion with center and surround kernels h2, and h f in Equation 60 
9), masking (via the masking kernel h m , and Equation 11), 
and non-linear pooling over space (via the power function and 
pooling convolution in Equation 12). Masking can provide an 
important function, because the detection of blur comprises 
detection of a small contrast (the departure from the perfect 65 
edge) superimposed on a large contrast pattern (the edge 
itself). 


symbol 

definition 

units 

example 

value 

usable 

range 

s c 

center scale 

degrees 

2.77/60 

±50% 

s s 

surround scale 

degrees 

21.6/60 

±50% 

s m 

masking scale 

degrees 

10/60 

1/60 - 60/60 

T 

masking threshold 

contrast 

0.3 

0-1 

S 

sensitivity 

dimensionless 

217.6 

±50% 

P 

pooling exponent 

dimensionless 

2 

1-6 

K 

adaptation weight 

dimensionless 

0.772 

0-1 


It will be understood that the descriptions of one or more 
embodiments of the present invention do not limit the various 
alternative, modified and equivalent embodiments which may 
be included within the spirit and scope of the present inven- 
tion as defined by the appended claims. Furthermore, in the 
detailed description above, numerous specific details are set 
forth to provide an understanding of various embodiments of 
the present invention. However, one or more embodiments of 
the present invention may be practiced without these specific 
details. In other instances, well known methods, procedures, 
and components have not been described in detail so as not to 
unnecessarily obscure aspects of the present embodiments. 
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What is claimed is: 

1 . A method of measuring visual motion blur *P comprising 5 
obtaining a moving edge temporal profile i 1 (k) of an image of 
a high-contrast moving edge with at least 10 samples across 
the step in r^k); 

calculating m 1 (k)=m(k) for r(k)=Ti(k), and 

m 2 (k)=m(k) for r(k)=r 2 (k), wherein r 2 (k) is the wave- 10 
form for an ideal step edge waveform with the same 
amplitude as r x (k); and 
calculating W as a difference function, 


visual motion blur is based on obtaining the moving 
edge temporal profile rl(k) of an image, calculating 
ml (k) and m2(k), and calculating visual motion blur as 
a difference function. 

2 . The method of claim 1 , wherein the center kernel is given 
by 


1 kLx 
h c {k) = — sech/r . 


3. The method of claim 2, wherein the surround kernel is 
given by 


where Ax is the sample interval in degrees of visual angle, 
S and (3 are parameters, 

m(k) is the masked local contrast, 



m(k) = 


c(k) 

V 1 + e(k) 


4. The method of claim 3 wherein the surround kernel scale 
s 5 is about 7.8 times the center kernel scale s c . 

5. The method of claim 2, wherein the masking kernel is 
given by 


e(k) is the effective local contrast energy, 


e(k) = h m (k)^ | 2 , 

30 

h,„(k) is a masking kernel having a masking scale s m and 
unit area, 

T is a masking threshold parameter, 
c(k) is the local contrast waveform, 


h c (k)*r(k ) 

c(k ) = — — — = - 1, 

K[h s (k) * r(k)] + (1 - k)R 

h c (k) is a center kernel having a center scale s c and unit 
area, 

U(k) is a surround kernel having a surround scale s 5 and 
unit area, 

R is the mean relative luminance of r^k) or, 
k is an adaptation weight parameter, and 
r(k) is an arbitrary waveform; 

wherein the offset of the ideal step edge relative to r^k) is 
adjusted to minimize tR; utilizing the measured visual 
motion blur to provide just noticeable differences of 
motion blur in display technology, where the measured 



6. The method of claim 5 wherein the masking kernel scale 
s„, is about 3.6 times the center kernel scale s c . 

7. The method of claim 5 wherein the center kernel scale s c 
is between about 1 .38/60 and about 4.1 6/60 degrees of visual 
angle, the surround kernel scale s i; is between about 21.6/60 
and about 32.4/60 degrees of visual angle, and the masking 

35 kernel scale s„, is between about 1/60 and 60/60 degrees of 
visual angle. 

8 . The method of claim 7 , wherein the center kernel scale s e 
is about 2.77/60 degrees of visual angle, the surround kernel 
scale s 5 is about 10.8/60 degrees of visual angle, and the 
masking kernel scale s m is about 10/60 degrees of visual 

40 angle. 

9. The method of claim 1, wherein T is between 0 and 1 . 

10 . The method of claim 9, wherein T is about 0.3. 

1 1 . The method of claim 1, wherein S is between about 1 09 
and about 326. 

45 12 . The method of claim 11 , wherein S is about 217.6. 

13 . The method of claim 1 , wherein |3 is between 1 and 6. 

14 . The method of claim 13 wherein [3 is 2. 

15 . The method of claim 1 , wherein k is between 0 and 1 . 

16 . The method of claim 15 , wherein k is about 0.772. 





